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For phase transitions in disordered systems, an exact theorem provides a bound on the finite 
size correlation length exponent: ups > 2/d. It is believed that the true critical exponent of a 
disorder induced phase transition satisfies the same bound. We argue that in disordered systems the 
standard averaging introduces a noise, and a corresponding new diverging length scale, characterized 
by vfs ~ 2/d. This length scale, however, is independent of the system's own correlation length ^. 
Therefore v can be less than 2/d. We illustrate these ideas on two exact examples, with f < 2/d. 
We propose a new method of disorder averaging, which achieves a remarkable noise reduction, and 
thus is able to capture the true exponents. 
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Using a very general formulation, Ref. presented 
an exact theorem, which puts constraints on the finite 
size correlation length exponent i^ps of a large class of 
disordered systems: vps ^ 2/(i, where d is the dimen- 
sion. This relation is often referred to as the quantum 
Harris criterion ||^ . While many investigations found ex- 
ponents in accordance with this bound, there is an in- 
creasing number of results in contradiction with it. In 
particular, in a model for charge density waves exact cal- 
culations yielded = 1/2 below four dimensions ||^, and 
numerical studies on 2d disordered Bose-Hubbard mod- 
els found v ~ 0.7 [Q. Experimentally the Bose glass 
transition of helium in aerogel ||^, and the localization 
transition in doped semiconductors seem to violate 
this bound. In this paper we argue that the vps > 2/d 
constraint is characteristic only to the method used to 
carry out the disorder average, and the true exponent v 
is independent of this bound. 

To start our considerations of random sytems, we chose 
the same type of disorder used by Ref. [|l| : a binary distri- 
bution for, say, a disordered site energy. Typically, phys- 
ical quantities are calculated by averaging over different 
disorder realizations. For calculational convenience, the 
standard method is analogous to the "grand canonical" 
approach: impurities are put on each site with a given 
probability, p, and the averaging is carried out for all 
possible concentration of impurities and their configura- 
tions. An alternative method, which could be termed 
the "canonical" approach, keeps the number of impuri- 
ties fixed, and the average is taken only over the possible 
configurations of these impurities. For infinite systems 
the two methods are equivalent. The density fluctua- 
tions in the grand canonical method, however, introduce 
an extra noise. This noise vanishes in the infinite sys- 
tem, but it may alter the results of the finite size scaling. 
The "canonical averaging" strongly reduces this noise by 
excluding density fluctuations. 

We now argue that the bound obtained in Ref. is 
only generated by the noise introduced by the "grand 



canonical averaging" . Different choices, such as using 
"canonical averaging" , produce different bounds. The 
theorem of Ref. Q considers a random system where 
a phase transition is induced by changing the concen- 
tration K of site (or bond) impurities. Let Y be any 
event depending on disorder realizations in a finite vol- 
ume, with probability V{K). This 'P{K) is calculated by 
averaging over all disordered configurations, and select- 
ing those compatible with Y. Averaging is performed 
in the "grand canonical" way, since fiuctuations in the 
density of impurities are allowed. From these premises 
the exact statement \dV{K)/dK\ < const. \/N follows, 
where N is the system size. A closer look at the proof 
reveals that this result is derived solely from the con- 
centration fiuctuations of the impurities, which were ex- 
ternally introduced in the averaging process (see the last 
equation of the proof in Ref. [Q). Thus the bound on 
\dV{K)/dK\ does not relate to the intrinsic properties 
of the system under investigation. It only reflects the 
"resolution" of the "grand-canonical averaging" . In other 
words, because of the presence of the density fluctuations, 
the minimal resolvable change in K is dK cx 1/ ^/N . The 
probability V can change at most 0(1), immediately ex- 
plaining the above bound. 

On the other hand, if one uses "canonical averaging" , 
then the above inequality does not apply. For, in con- 
trast to the previous case, the number of impurities is 
now well defined. In the present binary example, the re- 
solvable change of K is bounded only by its minimum 
allowed increment, 1/iV. Hence, \d'P{K)/dK\ < N . 
Along the lines of Ref. |]l|, the inequality vps > 1/rf now 
follows. As before, this inequality is characteristic of the 
"canonical averaging" only, and does not impose any re- 
striction on the true exponent u of the physical system. 
The physical reason behind this is that both averaging 
procedures introduce a new characteristic length scale, 
which has the potential to obscure the true correlation 
length of the physical system. 

It is also important to note that the assumption of a 
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binary disorder plays a crucial role in deriving the above 
bounds. For continuous distributions they do not nec- 
essarily apply. To see this, consider the following sim- 
ple example, motivated by the quantum phase transi- 
tion between the so called Mott-Insulator and Bose-Glass 
phases, which takes place in interacting bose systems 
with site disorder. At this transition the renormaliza- 
tion group flows are controlled by a fixed point with zero 
hopping strength , thus the system reduces to a collec- 
tion of independent sites with random energies. Let the 
distribution of the site energy e € [0, K] be 

whith a > — f. We generate N independent {i — 
1,...,N) from the above distribution. We define the 
finite-size event Y to occur, when all e^'s are smaller then 
a given value /i € {0,K]. We fix the value of fi, and 
drive the transition by changing K . As required by the 
theorem of |jl|] , the probability VofY happening is finite 
at the critical value of the disorder, Kc ~ /i. It goes to 
zero exponentially with the system size N for K > Kc- 
Close to the transition, for S ~ {K — Kc)/Kc ^ 1, this 
probability is 

ViN,S)^e-^^^^' . (2) 

A characteristic length scale ^/ can be now defined as 
a function of 6. It is determined from the system size 
as Nf = ^f, where V{Nf,5)/V{Nf,Q) - I/e. Defining 
a critical exponent as ^/ oc 5~'^''^ one arrives at vps = 
(a -|- For a < 1, vps is less than 2/d. While we 

considered a concrete example, we emphasize that this 
result can be relevant for any transition driven by local 
singularities in the action. 

Motivated by the above observations, we now attempt 
to construct a modified finite size scaling procedure, which 
does have the potential to access the true exponents of 
the system. The centerpiece of our argument is the fol- 
lowing observation. If the distribution of the disorder is 
given in an analytic form, that uniquely determines 
the critical value of the control parameter for the infinite 
system. However, any given disorder realization in a fi- 
nite system could have been generated from disorder dis- 
tributions with a range of parameters, corresponding to a 
range of Kc values. In other words, it is unclear which in- 
finite system's finite size realization did one simulate. A 
distribution is characterized completely by its moments. 
Typically, Kc is linked to some of these moments, for in- 
stance the dispersion. For a finite system of size TV, this 
dispersion is determined only with a relative uncertainty 
of 0{1/\/N). Therefore there is a range of distribution 
parameters which are compatible with the specific real- 
ization, and thus could have generated it. This raises the 
problem, which Kc to use in a finite size scaling analysis. 



The standard procedure answers this question by as- 
suming that one can use a single Kc for all samples gen- 
erated from the same distribution. However, the above 
argument suggests that the very same sample may be 
the realization of distributions with different parameters, 
leading to an inherent noise in the procedure, similar to 
the above considered binary examples. In order to avoid 
such a built-in noise, we now propose a modified finite 
size scaling procedure for disordered systems. We sug- 
gest that for each disorder realization one should identify 
the distribution, and in particular the critical parameter 
Kc, which it most likely corresponds to. In practice this 
might be difficult, and we return to this question later. 
For the moment, we only assume that it is possible to 
identify Kc- We propose that the natural control pa- 
rameter of the critical behaviour is A = {K — KD/K"^. 
The act of averaging should then be performed for the 
samples with the same A. We propose to adopt the fi- 
nite size scaling hypothesis for the critical behaviour of 
a generic physical quantity Q, 

QiL,A)^L-yqiLAn , (3) 

where q{z) is a universal scaling function, and y, v are the 
critical exponents for Q, and the true correlation length 
^ cx A^''. Note that some aspects of this proposition 
are already practiced in numerical studies: sizeable noise- 
reduction is customarily reached by adjusting the random 
variables after they are generated, e.g. in order to keep 
their mean value constant. 

Next we assume the validity of Eq.|| and perform the 
standard finite size scaling, to demonstrate how that pro- 
cedure's inherent noise can mask the true critical be- 
haviour. Some of the key results of the analysis are: i) we 
find that the exponent of the intrinsic correlation length 
v might be different from vps appearing in the standard 
finite size scaling. Therefore the theorem of Ref. ^ does 
not provide constraints on the true exponent v. ii) In 
particular, v can he less than 2/d. In this case typically 
VPS = 2/d. 

The standard finite size scaling procedure in dis- 
ordered systems calls for calculating a physical quantity, 
Q, such as the critical susceptibility, for different val- 
ues of N and K, the system size and control parameter, 
each time performing the calculations for a number of 
disorder realizations. Averaging over the disorder yields 
{Q{K)), and the critical coupling Kc is then identified 
for instance from a crossing pattern Requiring the 
collapse of the data, when plotted as a function of L^^'^S, 
where S — {K — Kc)/Kc, determines the exponents. 

To make contact between the standard scaling proce- 
dure and Eq.^, a relation between the unique Kc and the 
fiuctuating K"^ has to be constructed. A simple represen- 
tation of the inherent noise, or unccrtainity, is to assume 
the validity of the central limit theorem for Kc 
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where x is a random variable with a distribution width 
of 0(1). Here D measures the scatter in X^, and 5 is the 
distance from the average critical point Kc- As we have 
seen, this is not necessarily true for all systems, and we 
will return to the case when the fluctuations scale with 
a different power. 

The standard procedure neglects the fluctuations of 
if^, which is equivalent to averaging Q over the random 
variable x of Eq.|[ 

{Q) = L-yl^ql^D'^L^-^{x+^J^yy^ . (5) 

Here the x average is denoted by (...), corresponding to 
the standard averaging procedure, as opposed to Q, the 
correlated averaging of the new procedure in Eq.^. 

First we analyze the critical point itself, then we shall 
proceed to extract the critical behaviour of the correla- 
tion length. At (5 = the scaling form for Q is 

{Q)^L-y(q(D'^x'^L^-'^)) . (6) 

For V > 2/d the argument of the scaling function ap- 
proaches zero with increasing system size, and the L de- 
pendence of the averaged quantity {Q(L)) is characterized 
by the intrinsic exponent y. Here we use the customary 
assumption that the universal scaling function q[z) ap- 
proaches a finite value as z ^ 0. 

In the u < 2/d case, however, the argument of q{z) 
goes to large values, probing deeply non-critical regions, 
even though the system is assumed to be at criticality. 
To highlight the consequences of this, we proceed with a 
generic form for the asymptotic behaviour of the scaling 
function, adopting q{z) cx z' "'3. From Eq.| (Q) 
where 'y = y + /3(l — dv/2). Clearly the L dependence of 
the averaged (Q) is different from the intrinsic value y. 

Next we develop an understanding of the region in the 
proximity of the critical point, i.e. the case of finite S. Let 
us first focus on i/ < 2/d. From Eq.^ one identifies two 
scaling regions, governed by two different characteristic 
diverging length scales. 

For large system sizes inevitably D"^ L^^'^'^^^ 1 , so 
the argument of q{z) again extends to large values. Uti- 
lizing the previous asymptotic model form, 

{Q)^L-'^q{SL^/^) . (7) 

from which a length scale can be identified, characteriz- 
ing the finite size scaling of (Q), averaged in the stan- 
dard way. It diverges with an exponent vps = 2/d even 
though the true exponent v is less than 2/d. This result 
now demonstrates in general, what has been observed 
earlier for the binary example: the standard, or "grand 
canonical" averaging introduces a noise, which in turn 
generates a new length scale and a corresponding new 
exponent into the analysis. 



The other scaling region is reached when SL'^^^/D ^ 1 . 
In this limit 

{Q)=L-yq{6''L) . (8) 

As is known, for large values of S'^L, the exponent is 
not accessible by finite size scaling hence S'^L should 
be kept around unity. Therefore the determination of v 
requires the study of the region away from the asymp- 
totics: large 5 and small system sizes. For weak disorder 
{D <^ 1) this window in fact might be wide enough for 
practical purposes. To reiterate, however, studies con- 
centrating on the asymptotic region are bound to see 
VPS = 2/cf. 

In the case oi v > 2/d the standard procedure is capa- 
ble of accessing the true v. it can be obtained from (Q) by 
increasing the system size to the extent of 5L'^^^ / D ^ 1, 
but keeping cx 0(1). This again implies avoiding 
the "non-scaling" region around 5 = 0. For strong dis- 
order and small available system sizes, one can end up 
again with large arguments of q{z), and consequently in 
the scaling regime described by vps ^^nd 7 (Eq.^. There 
are several additional crossover regimes in the parameter 
space, which can be studied based on Eq.^. 

What happens if instead of the central-limit-theorem 
form, L^"^/^, the fluctuations olK"^ are described by some 
other power law? For instance, on physical grounds, the 
fluctuations may scale with the correlation length expo- 
nent 



Substituting this expression into Eq.|| and averaging over 
X shows that the standard procedure and our proposition 
give the same result for the exponents, although the scal- 
ing function changes due to the difference in averaging. 

Equation ^ can help incorporate the idea of in the 
scaling. One can appreciate that from looking at a finite 
sample, it is far from trivial to identify K^, belonging to 
the infinite system. A solution might be suggested by re- 
calling that for ordered classical magnets, the maximum 
of the susceptibility of a finite size sample is shifted as: 
Tc(L) — Tc(oo) (X L"^/"^, just as in Eq.^ Scaling then can 
be performed in terms of Tc(L), resulting in the correct 
exponents. Even in the absence of knowledge of the 
of the corresponding infinite system, one can still extract 
a from a specific feature of a critical quantity of the 
finite size system. Using this in our new scaling ap- 
proach should provide the correct exponent i^, provided 
that — cx L~^^'^ , a reasonable assumption. 

We are thus left with the task of identifying if^ of a 
finite system. For many quantum systems at T = a 
reasonable proposition for might be the value of K , 
where the gap to the first excitation vanishes or has a 
minimum. For classical systems may be identified 
where some critical susceptibility exhibits a maximum. 
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To demonstrate the above ideas, consider strongly in- 
teracting bosons in a random potential at zero temper- 
ature. In Ref. renormalization flows were generated 
by integrating out the sites with highest excitation en- 
ergies. For infinite range hopping the renormalization 
group (RG) equations are exact. In particular, at the 
Mott-Insulator to Superfluid transition weak disorder is 
irrelevant and i' — 1/d. Because of the presence of an 
underlying RG one expects the validity of finite-size scal- 
ing. 




FIG. 1. Scaling plot of the inverse susceptibility using the 
novel and the standard (insert) averaging procedure for sys- 
tem sizes N=64, 128,256. 

We carried out the finite-size scaling analysis of the 
average local susceptibility at weak but finite disorder for 
system sizes N = 64, 128, 256. First we used the standard 
averaging procedure (insert of Fig.l), and we obtained 
i/pg ~ 3/d after averaging over 1024 realizations of a 
uniform disorder distribution of the random potential. 
The collapse of the curves for different system sizes is 
not perfect, and we expect that as N increases, vps 
2/d. Fig.l shows the same quantity scaled by using KJ^ 
extracted from the divergence of the susceptibility for 
each sample separately. The scaling is convincing, and 
yields the exact exponent v — 1/d. The exhibited curves 
were obtained by averaging over much fewer samples than 
before, only 16, yet the scaling region extends by more 
than an order of magnitude further in terms of the scaling 
variable, NA, clearly demonstrating a very effective noise 
reduction. 

In some numerical studies, such as in Rcf. Q a vps < 
2/d has been reported, using the traditional averaging 
procedure. We would like to emphasize that this find- 
ing can be perfectly accomodated in the present theory. 
First, our analysis does not suggest that vps must be 
greater or equal to 2/d: this is only the most likely sce- 
nario. If, for instance, Eqjo] describes the fluctuations of 
K^, then i/ps = and thus can be less than 2/d. Appar- 
ently, this is the case in the example of the Mott-Insulator 
to Bose-Glass transition in Eq||. Second, as emphasized 



after Eq.|[ if the fluctuations of K'^ are small and the 
sample size is not too big, then the true v can and will 
be observed in finite size scaling. Finally, this theory 
is not addressing the problems associated with distribu- 
tions with long power-law tails pO[ | , or multicritical fixed 
points [0. After averaging, we still expect a gaussian 
distribution, with exponential tails. 

Now we would like to reflect on the Harris criterion . 
An insightful derivation imagines dividing the sample to 
blocks of size ^. For v <2/d, the fluctuations of the local 
"Tc's" of the blocks are bigger than the distance from the 
true Tc and it is concluded that u cannot be smaller than 
2/d. In our framework this negative result only means 
that the blocking procedure ceases to be a valid approach 
to the infinite system, and says nothing about the value 
of the inherent exponents. An RG based investigation of 
this problem will be given elsewhere. 

In sum, we reinvestigated the theory of finite size scal- 
ing in disordered systems. We found that the standard 
averaging procedure introduces a new diverging length 
scale into the problem, therefore the finite size scaling 
exponent vps may be unrelated to the inherent v of the 
true correlation length. In particular, we constructed 
two examples explicitly, where exact calculations proved 
that the inherent v is less than 2/d. We proposed an 
alternative method, which achieves a remarkable noise 
reduction, and therefore is capable of accessing the true 
exponents of the physical problem. 
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ported by NSF DMR-95-28535. 



J. T. Chayes, L. Chayes, D. S. Fisher, and T. Spencer, 

Phys. Rev. Lett. 57, 2999 (1986). 

A. B. Harris, J. Phys. C 7, 1671 (1974). 

O. Narayan and D.S. Fisher, Phys. Rev. B46, 11520 

(1992). 

J. Kisker and H. Rieger, |cond-mat/9612027 . 



P. Crowell, F. Van Keuls and J. Reppy, cond 



mat/9612133 



[6] 



[7] 
[8] 

[9] 
[10] 

[11] 



M. Paalanen, T. Rosenbaum, G. Thomas, and R. Bhatt, 
Phys. Rev. Lett. 48 1284 (1982); H. Stupp et al, Phys. 
Rev. Lett. 71, 2634 (1993). 
F. Pazmandi and G. T. Zimanyi, 



:ond-mat/9702169 



See the papers in Finite Size Scaling, ed. J.L. Cardy, 

North-Holland (1988). 

K. Binder, Z. Phys. B43, 119 (1981). 

D. Fisher, Phys. Rev. B 51, 6411 (1995). 

R. Singh and M. Fisher, Phys. Rev. Lett. 60, 548 (1988). 



4 



